A single ion as a three-body reaction center in an ultracold atomic gas 
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We report on three-body recombination of a single trapped Rb+ ion and two neutral Rb atoms in an ultracold 
atom cloud. We observe that the corresponding rate coefficient A'3 depends on collision energy and is about 
a factor of 1000 larger than for three colliding neutral Rb atoms. In the three-body recombination process 
large energies up to several 0.1 eV are released leading to an ejection of the ion from the atom cloud. It is 
sympathetically recooled back into the cloud via elastic binary collisions with cold atoms. Further, we find that 
the final ionic product of the three-body processes is again an atomic Rb+ ion suggesting that the ion merely 
acts as a catalyzer, possibly in the formation of deeply bound Rb2 molecules. 

PACS numbers: 34.50.-s, 34.50.Lf, 37.10.-x, 37.10.Rs 



Early on in the quest for ultracold quantum gases three- 
body recombination played a crucial role as a limiting fac- 
tor for Bose-Einstein condensation. It was first investigated 
in spin-polarized hydrogen |ll|] and somewhat later for alkali 
atoms Qlstl- For samples with large scattering lengths three- 
body processes can be resonantly modified through the forma- 
tion of Efimov states [4]. Recently, three-body recombination 
was investigated with single atom resolution ISJ. Combining 
ultracold atoms with cold trapped ions is an emerging field 
where large scattering cross sections naturally come into play 
due to the comparatively long range l/r"* polarization interac- 
tion potential. Two-body collisions between atoms and ions 
in the low energy regime have been recently studied ||6- 121. 
In this letter, we report on three-body collisions involving two 
ultracold **^Rb atoms and a ^^Rb+ ion at mK temperatures. 

The ion in our experiment can be regarded as a reaction 
center, facilitating molecule formation through its large in- 
teraction radius. For the work presented here, it is essential 
that we work with ions and atoms of the same species. This 
renders charge transfer reactions irrelevant, which otherwise 
would strongly constrain our measurements. Since no acces- 
sible optical transition is available for Rb+, it is not amenable 
to laser-cooling and cannot be imaged. We therefore detect 
the ion and investigate its dynamics in an indirect way, i.e. 
through its action on the atom cloud. In our experiments, we 
place a single ion into the center of an atomic sample result- 
ing in a continuous loss of atoms due to elastic atom-ion col- 
lisions. This behavior is interrupted when a highly energetic 
three-body process ejects the ion from the atom cloud. By 
examining the statistics of ion-induced atom loss in hundreds 
of repetitions of the experiment, we can investigate a num- 
ber of important details of the three-body process, such as its 
quadratic density-dependence, the energy that it releases, its 
rate coefficient Kt, and the dependence of Kt, on collisional 
energy. We also obtain information on the reaction products. 
As an important byproduct, our measurements also demon- 
strate sympathetic cooling of an ion from eV energies down 
to about 1 mK using an ultracold buffer gas. 

The atom-ion collision experiments are conducted in a hy- 
brid apparatus (for details see Ill3ll ) where a sing le **^Rb+ ion, 
trapped in a linear Paul trap, is brought in contact with an ul- 



tracold cloud of spin polarized ^^Rb atoms {F = l,m/r = — 1). 
The atom cloud is previously prepared at a separate location 
from where it is transported to the Paul trap and loaded into 
a far off -resonant crossed optical dipole trap. The dipole trap 
is at first spatially separated from the trapped ion by about 
50 /xm. To start the atom-ion collision experiments it is then 
centered on the ion with /im precision within a few 100 ms. 
At this point the atom cloud consists of A^at ~ 4.0 x 10^ atoms 
at a temperature of Tat « 1 .2 /xK and a peak density Mat « 
1.1 X 10^^cm"\ At trapping frequencies of (190, 198,55)Hz 
this results in a cigar shaped cloud with radial and axial exten- 
sions of about 10 ^m and 35 /xm, respectively. 

The single Rb+ ion is confined in a Paul trap driven at a 
frequency of 4.17 MHz resulting in radial and axial trapping 
frequencies of 350 kHz and 72 kHz, respectively. As the trap 
is about 4 e V deep, the ion remains trapped on timescales of 
days and can typically be used for thousands of experimental 
cycles. It is initially produced by photoionization of an atom 
from a cold Rb cloud in the Paul trap |14]. Typical kinetic 
energies Eion of the ion after sympathetic cooling in the atom 
cloud are about a few mK-^B where is the Boltzmann fac- 
tor This energy scale is mainly set by two quantities: (1) The 
excess micromotion (eMM) 1 15] in the Paul trap whose main 
part we can control by compensating stray electric fields lll4[l . 
(2) Heating effects induced by the interplay of micromotion 
and elastic collisions Il6i418ll . 

As described in fl, an ion immersed in an ultracold atom 
cloud leads to atom loss by expelling atoms from the shal- 
low optical trap (w 10/iK A;B trap depth) via elastic colli- 
sions. The radio frequency (rf) driven micromotion is a con- 
stant source of energy which drives these loss-afflicting colli- 
sions. Figure[T^ shows such a decay of an atom cloud at rela- 
tively low densities («10'' cm^^) and relatively high ion en- 
ergies (« 35mK-A:B ijj). Plotted is the number of remaining 
atoms after an atom-ion interaction time t. Each data point 
corresponds to a single measurement. Overall, the plot shows 
a relatively smooth decay of the atom cloud with a relative 
scatter of the atom number of less than 10%. This changes 
drastically when we carry out the experiments at low ion en- 
ergies (»i0.5 mK-^B ilStl) and larger densities (sslO'^cm^^) 
(Fig.llJ)). Here, the scatter dramatically increases with T and 
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FIG. 1 . Decay of the atom cloud under influence of a single trapped 
ion. (a) Remaining atom numbers after interaction time T for an 
ion with ~ 35mK-fcB 1 19] and n^t ~ lO" cm^-'. The solid line 
indicates the decay of the mean atom number, (b) Same as (a) but 
Eioa ~ 0.5mK A:B fisll and nat ~ lO'^cm^-'. (c) Histograms of the 
data shown in (b). The vertical axis counts the number of incidents 
for a particular atom number A'at and interaction time T. 



is on the order of the number of lost atoms. In Fig.[T]; his- 
tograms of the data in Fig.[lJ) are shown. With increasing time 
T the initial normal distribution becomes bimodal as a strik- 
ing tail towards large atom numbers emerges. At the tips of 
the tails we find cases where even after interaction times of 
several seconds barely any atoms have been lost, a signature 
of missing atom-ion interaction. Apparently, sporadically the 
ion is ejected from the atom cloud and promoted onto a large 
orbit for a period of time during which atom-ion collisions 
are negligible (Fig.|2^). In principle, this is reminiscent of 
the energy distributions with high energy tails that have re- 
cently been predicted for trapped ions immersed in a buffer 
gas [16, 17]. However, it turns out that such an explanation is 
inconsistent with our observations lEoll . Rather, we find that 
it is a three-body recombination process involving the ion and 
two neutrals that ejects the ion from the cloud. If, for example, 
a ground state molecule is formed, binding energies on the or- 
der of 0.5 eV can be released. Due to the large trap depth the 
ion is not lost in such an event, but it is recooled back into the 
cloud through binary collisions after some time. 

Figure [SJ) illustrates in a simple picture how the decay of 
the atom number over time can follow different paths. The 
solid trace Ti shows the case when only binary atom-ion col- 
lisions occur Such traces result in the narrow peak of the atom 
number distribution sketched on the right of Fig.|2i). Traces T2 
and T3 exhibit three-body collisions at points E2 and E3. At 
point R2 the ion reenters the atom cloud after recooling. Rare 
three-body events and long interruption times Tout result in a 
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FIG. 2. (a) Illustration of an atom-atom-ion collision. (I) Two 
atoms simultaneously enter the interaction radius of the ion and a 
three-body process takes place. (II) The three-body reaction ejects 
the ion onto a trajectory much larger than the atom cloud, (b) Il- 
lustration of our simple model. Left: Various possible time traces 
for the atom number. If only binary atom-ion collision occur the 
atomic sample decays exponentially (Trace Ti). Three-body events 
(E2, E3) interrupt the atom loss until the atom is recooled and reen- 
ters the sample at point R2 (Traces T2 and T3). Right: Sketch of 
the resulting atom number distribution when averaging over many 
experimental runs, (c) Plot of the probability /threebody for initial 
atomic densities (1.8,1.1,0.7,0.3) x lO'^cm^^ and atom numbers 
(6.5,4.0,2.8, 1.6) X 10"*, respectively. The solid lines are results of 
the numerical simulation. 



long tail of the distribution. We can reproduce the histograms 
in Fig.[T]with a simple Monte Carlo type simulation. We as- 
sume an initial normal distribution of the atom number which 
then decays exponentially with the binary atom-ion collision 
rate A'2nat- Here, K2~5x 10^^ cm^/s is a rate constant given 
by the product of the elastic cross section and the ion velocity. 
A three-body event, occurring at a rate K^nlf, interrupts this 
decay for a period fout. As the ion can only be recooled by 
the atomic sample, we assume the rate for reentry of the ion 
into the atom cloud to be proportional to the number of atoms 
1/ (fout) ~ Nat/ Coat with Com being a constant depending on 
the trap parameters. In order to fix the model parameters and 
rate constants it is convenient to carry out a data analysis as 
follows. We sum up the number of events in the tail of the 
distributions shown in Fig.[T]for each t and divide by the total 
number of measurements to obtain the probability Pthreebody 
that at least one three-body process takes place within time 
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T. These data are plotted as the filled black circles in Fig.|2j; 
(data set with Wat ~ 1-1 x 10^^ cm^^). The corresponding er- 
rorbars mainly reflect the difficulty to unambiguously assign 
data points to the peak or the tail of the distribution where both 
parts of the bimodal distribution strongly overlap. Besides the 
data of Fig.[T];, Fig.|2}; also contains data at three additional 
atomic densities. All four data sets have in common that the 
number of three-body events first rapidly increases and sub- 
sequently levels off. The leveling off is due to the fact that 
the probability for a three-body reaction is strongly density- 
dependent. Closer inspection of the data sets reveals a pecu- 
liar feature. In the beginning of the interaction (T < 1 s) only 
very few three-body events are detected for the lower density 
samples. We explain this delay by an initial phase of sym- 
pathetic cooling of the Rb+ ion which experiences significant 
heating during the preparation time of the atom cloud (-30 s). 
In fact, we observe ion heating specifically during evaporative 
cooling where the rf is ramped from 70 to 3 MHz. From nu- 
merical calculations we estimate that recooling times of about 
1 s in atom clouds with «at ~ 10^^ cm^^ correspond roughly 
to ion kinetic energies of a few 100K A;b. According to the 
calculations the ion will typically undergo several thousand 
binary collisions with cold atoms until it is sympathetically 
recooled to mK-^B energies. We are able to describe all four 
data sets in Fig.|2]; consistently with our simple Monte Carlo 
model (continuous lines). From independent fits to each data 
set we obtain nearly identical rate coefficients in the range 
= 3.1 — 3.5 X 10^^^ cm^/s. For each fit the initial cool- 
down time is accounted for by adjusting the starting time of 
the simulation. We note that the value for our atom-atom- 
ion rate coefficient is more than three orders of magnitude 
larger than the three-body coefficient for three colliding neu- 
tral ^^Rb atoms H . 

In order to challenge our analysis we have attempted to 
model the events that send the ion into orbit as two-body pro- 
cesses. The corresponding linear density dependence of the 
event rate yields much less consistent fit results. The two- 
body rate coefficients would differ by more than a factor of 3 
when comparing the analysis of the data set for the lowest and 
the highest density. As a cautionary note, we point out that 
three-body recombination processes to weakly -bound molec- 
ular states with binding energies < lOmeV are not detected 
in our experiments as the ion will not leave the atom cloud. 
Thus, the true three-body coefficient may even be significantly 
larger From our analysis we also extract an approximate 
value for the interruption time coefficient Cout ~ 1.7 X lO^s. 
For the typical atom numbers used here this results in sev- 
eral seconds of negligible atom-ion interaction following each 
ejection of the ion. 

In a further experiment, we quantify the kinetic energy 
AEion gained by the ion in a three-body event. The idea is to 
lower the depth of the ion trap such that an ion with an energy 
of a few 0.1 eV escapes while a cold ion remains trapped. The 
trap depth is reduced to one of 5 values J/red by lowering one 
of the endcap voltages of the Paul trap within 300ms. The 
voltage is held at this value for 200 ms and ramped back up 
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FIG. 3. Measuring the ion energy after a three-body process (black 
squares) and after purely binary collisions (grey dots). The trap depth 
is reduced to U^ed- (a) If the ion "survives" this procedure in the 
trap it will subsequently induce atom loss in a freshly prepared cloud 
of atoms. Plotted is the remaining atom number. Data points with 
^at lO"* signal the presence of an ion, while data points with A'at ~ 
4.5 X 10^ correspond to a lost ion. For better visibility, we have 
slightly offset in energy the black squares from the grey dots, (b) Ion 
loss probability Pjoss calculated from the data in (a). The continuous 
lines are guides to the eye. 



within 200ms. The trap depths f/red are indicated in Fig . 3. 
They have been computed using methods detailed in 112 III for 
both Rb+ (bottom abscissa scale) and Rbj (top). A negative 
trap depth value corresponds to non-trapping. After the trap 
depth reduction procedure is completed, the presence of the 
ion in the trap is probed via the loss it induces in a freshly 



prepared atom cloud containing about 5 x 10 atoms. For 
this probing procedure we deliberately apply an offset electric 
field of about 6 V/m to increase the eMM energy. In this way, 
we make three-body reactions unlikely and induce a rapid loss 
of atoms through binary atom-ion collisions. Figure |3^ shows 
the remaining atom number after 6 s of interaction time. An 
atom number < 1 x 10"* indicates the presence of an ion while 
a number around 4.5 x 10"*^ shows its absence. The clear split- 
ting of the two groups of data allows for ion detection with an 
efficiency close to unity. In addition, this ion detection is en- 
ergy resolved. Figure |3^ contains two different plot symbols, 
distinguishing two classes of ions that have undergone a dif- 
ferent prehistory before their energy resolved detection. Black 
plot symbols correspond to ions that have been promoted to a 
high energy orbit due to a three-body recombination. The re- 
combination occurred during 4s of atom-ion interaction and is 
detected through strongly reduced atom loss. Grey plot sym- 
bols correspond to ions where no suppression of atom loss 
was detected. These ions should in general have low kinetic 
energy. We now analyze the data points of Fig.|3^ by calculat- 
ing the probability for ion loss fioss (see Fig.[3]3). As expected, 
ions that were involved in a three-body recombination process 
can in general escape from deeper traps than cold ions. From 
the energy offset between the black and grey data we estimate 
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FIG. 4. (a): /fhreebody a function of the external electric field, (b): 
Double-logarithmic plot of /thieebody as a function of the ion energy 
^ion 1 19]. A scale for the three-body coefficients ^"3 as derived from 
the simulation is also given (see text for details). 



the gained energy A/Sion ~ 0.4eV. We note that for trap depths 
Ured < 0.25 eV the probabiHty of loss is high in general. This 
suggests that the stability of our trap is compromised at such 
shallow trapping potentials, limiting the accuracy with which 
we can determine the energy released in the three-body pro- 
cess. Still, we find a clear splitting between the black and grey 
data sets. Thus, a resolution of the measurement on the order 
of 0.1 eV seems plausible. 

Mainly two recombination processes come into consider- 
ation. In a reaction of the type Rb H- Rb H- Rb+ — > Rb2 + 
Rb+ the formation of a neutral molecule is catalyzed by the 
ion which carries away 2/3 of the energy released. If deeply 
bound Rb2 molecules are produced binding energies of up to 
~ 0.5 eV are released, in agreement with the measurement. A 
second possible recombination process, Rb H- Rb H- Rb+ — ^ 
Rbj + Rb, produces a molecular ion and a neutral atom. How- 
ever, as indicated in figure 3, the molecular ion, due to its 
higher mass, experiences a significantly shallower trap than 
the atomic ion and would immediately be lost for our parame- 
ter range. We thus infer that the ion at hand is Rb+. However, 
we cannot completely exclude the formation of an intermedi- 
ate molecular ionic state which may subsequently dissociate. 

In a third type of measurement we study the dependence of 
the three-body coefficient on the ion kinetic energy which we 
can tune by controlling the ion micromotion. For this we ap- 
ply a static electric field e perpendicular to the axis of the Paul 
trap and let the ion interact for T = 8 s with an atom cloud with 
«at ~ 1.0 X 10^^ cm^^. We find -Pthreebody to increase roughly 
by a factor of 5 when reducing e from 3.25 V/m to OV/m 
(Fig.|li). 

In order to express the electric field values in terms of ki- 
netic energy we make use of the relation Esmm — Ctmp ■ £^ + 
fires with Ctrap a Constant that depends on the trap configura- 
tion and the ion mass Iil5il . fires stands for residual uncom- 



pensated micromotion energy, e.g. due to rf phase delay be- 
tween the trap electrodes. The ion energy can be expressed as 
Eion = Cdyn ' £^eMM lfl9ll - Cdyn is a Constant which depends on 
the atom-ion mass ratio and the spatial extension of the atom 
cloud and for our experiments can be estimated to be about 2 
iITtIi . Since three-body coefficients often follow simple scal- 
ing laws, we attempt to describe our data with a power-law 
dependence of the form Kj, 0= E^^^ within our simulation. We 
have taken care to also account for the energy dependence of 
the two-body interactions. Good agreement with the data is 
achieved for a — —0.43, fires — 370/xK-^b and a maximal 
value for of 2.75 x lO^^^cm^/s (solid trace inFig.|4j3). For 
comparison, curves for exponents a = —0.5 and a = —0.33 
(dashed and dotted traces, respectively) are shown as well. 
A residual energy fires — 370/iK-^B is a reasonable value for 
our current setup and in agreement with other measurements 
of ours 114 1 . 

In conclusion, we have studied three-body recombination 
involving a single trapped ion and two of its parent atoms at 
collision energies approaching the sub-mK regime. With a 
relatively simple model we can understand the two- and three - 
body collision dynamics and extract corresponding rate coef- 
ficients. We observe an increase of the three-body rate coeffi- 
cient with decreasing collision energy, a behavior that can be 
expected to become crucial for future experiments targeting 
even lower temperatures. After a three-body event, ion en- 
ergies on the order of 0.4 eV were measured, indicating that 
deeply bound molecules have been created. Since we have 
not observed Rb^ ions, the formation of Rb2 seems probable. 
The ion would then act as an atomic size catalyzer at mK tem- 
perature. It is an interesting question whether such catalyzing 
action could also be observable close to more massive objects 
such as carbon nanotubes |22|]. Finally, as a byproduct of our 
investigations, we also observe sympathetic cooling of ions 
from energies in the 0.1 eV range back to mK temperatures. 
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